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Abstract. Muon pair production through the radiative return is of importance for a measurement of the 
hadronic production cross section in two ways: it provides an independent calibration and it may give rise to 
an important background for a measurement of the pion form factor. With this motivation the Monte Carlo 
event generator PHOKHARA is extended to include next-to-leading order radiative corrections to the 
reaction e + e~ — > fi + fJ,~j. Furthermore, virtual ISR corrections to FSR from pions are introduced, which 
extends the applicability of the generator into a new kinematical regime. Finally, the effect of photon 
vacuum polarization is introduced into this new version of the generator. 
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1 Introduction 



The high luminosity of currently operating and B- 
meson factories offers the unique opportunity to measure 
the famous ratio R = a(e + e~ — > hadrons)/<7 po i n t over a 
wide range of energies in one single experiment, through 
the radiative return. This possibility was suggested a long 
time ago for pion-pair production £Q . In view of the capa- 
bilities of the currently operating machines, the idea was 
revived in [2], where a Monte Carlo event generator EVA 
was constructed, which simulates the production of a pion 
pair plus one hard photon from initial- and final-state ra- 
diation (ISR and FSR). The program includes additional 
collinear radiation from the incoming electrons through 
structure function techniques and was used to demon- 
strate the feasibility of such a measurement in ongoing 
experiments. 



* Work supported in part by BMBF under grant number 
05HT9VKB0, EC 5th Framework Programme under contract 
HPRN-CT-2002-00311 (EURIDICE network), TARI project 
HPRI-CT-1999-00088, Polish State Committee for Scientific 
Research (KBN) under contract 2 P03B 017 24, BFM2002- 
00568, Generalitat Valenciana under grant GRUPOS03/013, 
and MCyT under grant FPA-2001-3031. 

a e-mail: czyzOus . edu . pi 

b e-mail: agrzel@us . edu . pi 

c e-mail: j ohann . kuehnOuni-karlsruhe . de 

d e-mail: german. rodrigoOcern. ch 



The generator was subsequently extended to simulate 
the production of other hadronic final states: four pions 
in and nucleon pairs in as well as muons 

The influence of collinear lcpton pair radiation on these 
measurements was investigated in [7]. 

An important step forward was the evaluation of the 
complete one-loop corrections to ISR [HIE] and the con- 
struction of the event generator PHOKHARA 0, which 
included these virtual corrections and the corresponding 
radiation of two photons. PHOKHARA, like EVA, can be 
easily used to simulate reactions with arbitrary hadronic 
final states, once a specific model has been chosen for the 
relevant matrix element of the hadronic current; it is cur- 
rently available for 27r, 47r, pp and nn production 
and of course for muon pairs. A number of cross section 
measurements based on the radiative return, which have 
made use of simulations based on EVA and PHOKHARA, 
have been presented recently: a precise measurement of 
the pion form factor made by the KLOE Collaboration 10 
I11II12II13) : preliminary results on four-prong final states 
have been obtained by the BABAR Collaboration |14j . 
Extensive discussions of theoretical and experimental as- 
pects of the radiative return, including various topics not 
covered by the present paper, can be found in the liter- 
ature |15Ul^lli7lll^lll0ll^0lliiill^ij|lijiilli>4ll^5llij^>lll>7llijaili>0ll^0| : 
considerations related to the radiative return can be found 
in (SUES; perspectives of its use outside hadronic cross 
section measurements and the long term potential of the 
method are outlined in 1331 . 
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Considering measurements with a precision close to 
one percent, FSR starts to become an issue, in particu- 
lar the contribution from the two-step process e + e~ — > 
77* (— > 7T + 7r - 7) JHj- FSR has to be modeled separately 
for every different mode. In leading order (LO), it is now 
included in PHOKHARA for the 7r+7r~ and the fi + fi~ 
modes, in next-to-leading order (NLO), i.e. for the two- 
step process with one hard photon from ISR and real plus 
virtual radiation from the final state for tt + it~ only. 

It is the purpose of the present paper to extend the 
same approach to final states with muon pairs. This is 
motivated, on the one hand, by the fact that this reac- 
tion is an important background for the pion form factor 
measurement, on the other hand, that it may be used for 
a measurement of the ratio between hadronic final states 
and muon pairs of the same invariant mass. This might 
lead to a direct determination of the i?-ratio, with the 
cancellation of many uncertainties that arise in the abso- 
lute cross section determination. 

Although, technically speaking, the treatment of muons 
is quite similar to that of pions, a number of physics as- 
pects are quite different. First of all, in contrast to radia- 
tion from muons, the amplitude for photon radiation from 
pions can only be parametrized by a model form factor, 
and the model dependence becomes particularly relevant 
for energetic photons. Then, at larger cms energies, rele- 
vant e.g. to the B-factories, the leading order FSR with 
only one photon in the final state is strongly suppressed 
in the pion case, as a consequence of the pion form factor, 
whereas the point-like behaviour of the muon leads to im- 
portant contributions from ISR as well as FSR in leading 
order. 

In order to arrive at a full result for muon pair produc- 
tion at NLO, corrections with combined virtual emission 
from the initial state and real hard emission from the final 
state (Fig. lc) must be included, and combined with those 
from "simultaneous" soft real ISR and hard FSR (Fig. la). 
To offer a uniform program, the same contributions will 
also be included for pion pair production. However, we will 
demonstrate that their effect is small already for KLOE 
energies, in particular for the standard KLOE cuts, and 
completely irrelevant for higher energies. 

Finally we study the influence of the leptonic and ha- 
dronic vacuum polarization on the radiative return, where 
we adopt the parametrization suggested in |34| . For the 
measurement of the integrated cross section with cuts sup- 
pressing FSR, vacuum polarization just leads to a multi- 
plicative factor that can be taken into account by correct- 
ing the result at the end of the experimental analysis. The 
forward-backward or the charge asymmetry, however, is 
affected by vacuum polarization effects in a non-trivial 
manner. 

Let us briefly outline the content of this paper. The 
amplitudes for real and virtual radiation in NLO for muon 
pair production, which are the main additional ingredients 
in the new program, will be introduced in Section 2. Tests 
of the technical precision of the program and some of its 
physics results will be presented in Section 3. The ISR cor- 
rection to real hard FSR for the 7r + 7r _ mode will be intro- 




(b) (c) 

Fig. 1. NLO corrections to the reaction e + e~ — > fi~"f: (a) 
double real emission, (b) virtual corrections to final-state ver- 
tex, and (c) virtual corrections to initial-state vertex. Only 
representative Feynman diagrams are depicted. 



duced in Section 4, which also contains a discussion of the 
influence of this additional term on present experimental 
studies. The implementation of vacuum polarization and 
its effect on the cross section and asymmetries will be the 
subject of Section 5. Section 6 contains a brief summary 
and our conclusions. 



2 FSR for muons at next-to-leading order and 
its implementation in PHOKHARA 

The most relevant contributions to FSR at NLO, which 
are depicted schematically in Fig.ri] consist of double real 
photon emission diagrams, where one photon is emitted 
off the initial-state leptons and the other is emitted from 
the final state (Fig. QJi), final-state vertex corrections to 
single initial-state photon emission (Fig.r^>), and initial- 
state vertex corrections to single final-state photon emis- 
sion (Fig.nk). 

The amplitude for the process depicted in Fig. 

e + (pi,A e +) + e~(p 2 ,K-) -> 

H + (qi,*»+) + (q 2 , X^-) + 7(fci,Ai) + 7(^2, A 2 ) , 

(1) 

is given by 

H (47ra) 3 / 2 / 

-^IFSNLo(A e + 5 A e - , A M +, A M -,Al, A 2 ) = | 

v\(pi,\ e +) A(A M +,A M -, Ai, A 2 ) u/(p 2 ,A e -) 

+ v\ I (p 1 ,X e +) S(A M +, A^-jAi, A 2 ) M//(p2,A e -)| 

+ (fo «-> fc 2 ,Ai «-> A 2 ) , (2) 
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where the index H indicates that both photons are hard, 
^(A M +,A M -, Ai, A 2 ) = 

(e*(fci, Ai)~fc+ - 2e*(k 1 ,X 1 ) ■ gi) £>-(A 2 , A M+ , A M - ) 
2fa ■ pi 

D-(A 2 , A M+ , V) (2e*(fci, Ai) ■ p 2 - fc+e*(fci, Ai)~) 
2fci • p 2 



+ 



(3) 



and 



B{\^+ , Xfj,- , Ai, A 2 ) — 

(e*(fci, AQ+fcf - 2e*(fci, Ai) ■ Pl ) £> + (A 2 , A M + , A M - ) 
2fa ■ pi 

£>+(A 2 ,A M+ ,A M -) (2e*(fci,Ai) -p 2 - fcf e*(fci, AQ+) 
2fci • p 2 



(4) 



with 



only, and 
77 V +V,£ 2 cut ) 



J+/3 'log(t) + l 



1 + /% 

0M 



-s' - 7ml + 3m^ar M 



log(2u/) 
- 2 + log 



1-/32 



2Li 2 



2/3„ 



1 + /3 M 



y-log(t)log , ( 



4 

1+^ 



(10) 



£" 2 ut is the maximal energy of the soft photon in the s' 
rest frame, and 



^ = .^1-47^2/5', t =- 



1-/3/, 



2ml 



(11) 



The virtual (Fig. [TJ;) plus soft photon corrections to the 
initial-state vertex can be written as |3f>| 



Q = pi + p 2 - fa = qi + q 2 + k 2 , s = Q . (5) 

D M is the current describing the final state 0], 

which we report below for completeness 



D" (A 2 , A M+ > A^-) = ie< ttj (g 2 , A M - ) A" (A 2 )u/ (<7i , A^+ ) 



+ A M - )B^(X 2 )v II {q 1 , A M +) 



(0) 



with 



A"(A 2 



(2 g2 ■ e*(fc 2 , A 2 ) + e*(fc 2 , A 2 )-fc+) o»- 
2k 2 ■ q 2 

a*- {2 qi -e*(k 2 ,X 2 ) + k+e*{k 2 ,X 2 )-) 
2k 2 ■ qi 



, (7) 



S^(A 2 ) 



(2 g2 ■ e*(fc 2 , A 2 ) + e*(k 2 , A 2 )+fc~) 
2k 2 ■ q 2 

qn+ (2gi ■ £*(fc 2 , A 2 ) + fc 2 -£*(fc 2 , A 2 )+) 
2fc 2 • qi 



, (8) 



and cr M± = (7,±(Tj), where Cj, i = 1,2,3, are the Pauli 
matrices and a — a^a^ for any four- vector a^. 

The virtual (Fig. plus soft photon corrections to 
the final-state vertex can be written as |35| 



^NLO,f = f »? V+ V> E^) daig( S >) , (9) 



where dojgj^ is the leading order e + e — * 7 cross 

section, with the photon emitted off the initial leptons 



acr IFSNLO, 



l = ^S v+s {8 ) ES ta )da^ R ( a ), (12) 



where dopsR is the leading order e + e~ — > /i + /j~7 cross 
section, with the photon emitted off the final leptons only, 
and, for ml <C s: 



S Y+S = 2|(L - 1) log (2w) + -L - 1 + C(2)| . (13) 

£"" lln is the maximal energy of the soft photon (or the 
minimal energy of the hard photon) in the s rest frame, 
while 



L = log — and w = E" 



(14) 



To match hard, soft and virtual radiation smoothly, the 
energy cutoff (i? 2 ut ) has to be transformed from the rest 
frame of the /i + /x~7 system (with the photon emitted from 
the final state) to the laboratory frame (e + e~ cms frame) 
(-E™ m ). In fact it is necessary to recalculate the soft pho- 
ton contribution, as the cut on Q 2 = (qi + q 2 ) 2 depends 
in the latter case on the angle between the two emitted 
photons. Now 



1 + 



log(2u>) + 1 + 

s' — s 



- 2 + log 



1-/3; 



log(*) 



~s' - 7ml + 3m^ p 



l + /3j 

/3u 



2Li 2 



2/3, 



£-]og(t)]og(i±& 



1 + /3 M 
(15) 



Using Eqs. J2J to (|15|l the implementation of FSR in com- 
bination with ISR is straightforward. 
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Table 1. Total cross section (nb), corresponding to the Q 2 
distributions from Fig. [5] for the process e + e~ — + n + /j,~'y at 
NLO, for different values of the soft photon cutoff. 

w y/s = 1-02 GeV 10.52 GeV 



10" 

10" 



28.237(1) 
28.236(1) 



0.19854(7) 
0.19883(14) 



We have not included diagrams where two photons 
are emitted from the final state, neither final-state ver- 
tex corrections with associated real radiation from the fi- 
nal state [37]. This constitutes a radiative correction to 
FSR and will give non negligible contributions only for 
those cases, where at least one photon is collinear with 
one of the muons. Box diagrams with associated real ra- 
diation from the initial- or the final-state leptons, as well 
as pentagon diagrams, are also neglected. As long as one 
considers charge symmetric observables only, their con- 
tribution is neither divergent in the soft nor the collinear 
limit and thus of order a/ir without any enhancement fac- 
tor. We want to stress that we have included only C-even 
gauge invariant sets of diagrams (see Fig. 9 of Ref. [T%| 
for a graphical representation of the equivalent set of di- 
agrams in pion pair production). The box and pentagon 
diagrams that we have neglected are related to ISR-FSR 
interferences, and therefore are suppressed after suitable 
kinematical cuts used within the radiative return method. 
One could also neglect some small contributions from the 
diagrams taken into account, but they are kept for the 
sake of completeness, as the full 1-loop radiative correc- 
tions will be consider in a future publication. At this point 
these small contributions will become relevant allowing ac- 
curate calculations not restricted to radiative return phys- 
ical configurations. 



3 Tests of the Monte Carlo program and 
discussion of physical results 

A number of tests were performed to ensure the tech- 
nical precision of the new version of PHOKHARA. The 
square of the matrix element, summed over polarizations 
of the final particles and averaged over polarizations of 
the initial particles, was calculated with FORM |2Hj, us- 
ing the standard trace method. External gauge invariance 
was checked analytically when using the trace method and 
numerically for the amplitude calculated with the helicity 
amplitude method. The two results for the square of the 
matrix element summed over polarizations were compared 
numerically. The code based on the result from the trace 
method was written in quadruple precision to reduce can- 
cellations. The code based on the result obtained with the 
helicity amplitude method uses double precision for real 
and complex numbers and is now incorporated in the code 
of PHOKHARA 4.0. Agreement of 13 significant digits (or 
better) was found between the two codes. The sensitivity 
of the integrated cross section to the choice of the cutoff 
w can be deduced from Fig. [21 and Tabled For simplicity 
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Fig. 2. Comparison of the Q 2 distribution for two values of 
the soft photon cutoff (w = 10" 4 vs. 10" 5 ) for y/s = 1.02 GeV 
(a) and yfs = 10.52 GeV (b). One of the photons was required 
to have energy > 10 MeV (for yfs = 1.02 GeV) and > 100 MeV 
(for y/s = 10.52 GeV). One sigma statistical errors are shown 
for the ratio of the two distributions. 



the same separation parameter w was chosen for ISR and 
FSR corrections. Choosing w — 10~ 4 or less, the result 
becomes independent of w, up to expected systematic dif- 
ferences of order 1CP 4 . These systematic differences are 
visible already in Fig. 2a, where x 2 /d.o.f = 30/14, while 
in Fig. 2b statistical errors are about 0.1%, and % 2 /d.o.f 
= 3.5/7. The tests prove that the analytical formula de- 
scribing soft photon emission as well as the Monte Carlo 
integration in the soft photon region are well implemented 
in the program. 

For the purpose of further tests we have also calcu- 
lated the integrated contribution of hard photon emission 
from the final state to the ISR spectrum. Integrating over 
all angles and energies, from £^ ut (defined in the s' rest 
frame) to the kinematical limit of the final-state photon, 
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we find that this contribution can be written as 

^i H FSNLO.f = - V H (s', £ 2 cut ) da[g ( S ') , (16) 

7T 



where 

+ C(2) 



1 + 



2 r 



Li 2 1 



- log(i m - t) 



U 2 (t m t) 



1+/9, 



Iog(*m) f 1 _ 



ft* 1 1 - $ 



(1-/%) 



1 



log 



17-/31- 2x » 



(1 - /9*)(1 - /&) 



log(* m ) 



;i J 33 - 17$ 



1 



1 - 3 2 



4/3 m /? M 
6(1-/32) 



-27+11/32-4-/ 



-21og(l-i m t) , 
i, and /J^ are defined in Eq. I|llfl . and 



(17) 



/3„ 



4m?, 



t, 



l-/3„ 
l + /3„ 



Here Q 2 ^ is the maximum value of Q 2 

Q 2 m = s ' - 2£ 2 cut V7 . 

For small -E™', ^ = E^/Vs' <C 1, the function ^ H 
reduces to 



,>',££"•) a log(2i») 



1 + ii Mo g(i ) 



A. 



+ 



log(*) 

8/3, 



Li 2 (t 2 )-C(2)-21og(i)log 



1 + /3 M 



-5 + /5 2 + 6x M 



2 log 



4$ 



11 

T 



(18) 



Adding virtual, soft (Eq. JTJJl) and hard (Eq. O) cor- 
rections, the familiar correction factor 39 



T]{S ) = 



+ 



log(i) log 
2 

/3^(3-/3 2 ) 



4Li 2 (i) + 2Li 2 (-t) 



(1+/3m) S 
8/3 2 

3fl, 



3 log 



1-$ 

4/3„ 



log(/3 p 



^(5-3/3 2 ) 31 ° g W ^o,oo,2 



48 



(33 + 22/3 2 - 7/3*) 
(19) 



is recovered. 
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Fig. 3. Comparison between the LO FSR cross section calcu- 
lated analytically (Eq. JSHJ) and calculated by PHOKHARA 
for a fixed value of w — E™ ln /^/~s. Errors represent one sigma 
statistical errors. 



For single-photon emission from the final state and no 
further photon radiation, a formula similar to Eq. <|16|) 
holds: 



FSR 



77 H ( S ,£f")a 



» > 



(20) 



with ?7 H defined in Eq. {TTJ) , £™ in defined now in the s 



rest frame and a 



(») 



the lowest order cross section 



e+e — — fj,~ 

of the process e + e~ — > 

The results of the tests are collected in Figs. |21 and 0] 
proving the excellent technical precision of this part of 
PHOKHARA, well below 0.05%. Results of similar tests at 
NLO are collected in Fig.0 where results of PHOKHARA 
are compared with the analytical result of Eq. (flfjf) . The 
latter checks the technical precision of the implementation 
of double-photon emission, where one photon is emitted 
from the initial state and the other from the final state, 
together with the corresponding virtual and soft correc- 
tions to the final-state vertex. An agreement of 1% on the 
function r\ means in fact an agreement better than 10 -4 
on the cross section thanks to the additional factor a/ir 
(compare Eq. © and Eq. J20J)- 

The relative size the NLO FSR contributions to the 
e + e~ — > /i + /i~7(7) cross section introduced in PHOKHA- 
RA depends on the event selection used. In the follow- 
ing, we will indicate some of its characteristic features. In 
Fig. the Q 2 differential cross section is shown with two 
peaks at low ('soft' muon pair production) and large Q 2 
(soft photon emission) values. In Fig. [(Jj the correspond- 
ing relative NLO FSR contributions are shown. They are 
<5 2 -dependent and might be as big as a few per cent. As 
seen by comparison of Figs. 03 and [(J: the relative size of 
this contribution does depend on the event selection used. 



At KLOE, the process e 



ground to the measurement of the 



7(7) is a back- 



cross 
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0.02 
0.015 

0.01 
0.005 


-0.005 

-0.01 h 
-0.015 

-0.02 



A 4 " 1 " A 4 7(7) 



ij( PHOKHARA) 


_ 1 






r\(ana 


lyiical) 
























♦ \ 

• 
















1 


























= 1.02 


GeV 

















0.2 



0.8 



0.4 0.6 

s'(GeV 2 ) 

Fig. 5. Comparison between analytical function 77 (Eq. 11910 
and T) obtained from PHOKHARA for = 1.02 GeV. 



emitted from the final state, are expected to be smaller, 
but will be relevant anyhow if the precision requirement 
for the background estimate is of the order of a few per 
cent. 

The process e + e~ — ► /i + 7 (7) is even more impor- 
tant for BABAR, where it is used as Q 2 -dependent lumi- 
nosity monitoring. As a result the accuracy requirement is 
higher. In Fig. [8] some features of the newly implemented 
corrections are shown in the interesting Q 2 range, from 
the point of view of using the radiative return method 
for measurements of the hadronic cross section. In Fig. [SJi 
the Q 2 differential cross section is given, including and ex- 
cluding NLO FSR corrections. As seen from Fig. |SJd the 
relative size of those corrections is quite small, of the or- 
der of 1%, if no angular cuts are applied; but, as shown 
in Fig. |Ht, their relevance depends on the event selection 
and the corrections can be as big as 10% for some event 
selection and Q 2 range. 



section by the radiative return method, due to possible 
pion muon misidentification. The angular cuts used by 
now by KLOE diminish the low Q 2 part of the cross sec- 
tion, while the track mass cut, which is specifically intro- 
duced to constrain the final state to configurations with 
7r + 7r~ and one photon only (see Ref. 0HI), does not allow 
for events with and one photon only (see Fig.U^, for 

corresponding cross sections) . If only angular cuts are ap- 
plied, the newly introduced corrections are at the level of 
a few per cent (Fig.03). However, if in addition the track 
mass cut is applied, leaving only events with two hard 
photons, the contribution of one-photon ISR plus one- 
photon FSR to the two-photon sample is up to 35% of the 
two-photon ISR sample (Fig. 0:), thus this contribution 
is indispensable for a reliable background estimate. The 
left-over two-hard-photon corrections, with both photons 



4 ISR corrections to FSR real hard emission 
for two-pion final state 

Initial-state vertex corrections to FSR for the reaction 
e + e _ — > 7r + 7r~7 (Fig. [30 were not included in version 
3.0 of PHOKHARA, and therefore, for consistency, dou- 
ble photon radiation (Fig. ^p) was restricted to emission 
of hard ISR. In the present version of the event generator, 
PHOKHARA 4.0, these corrections jSE] have been added 
and the generation of double photon emission is uncon- 
strained, provided that one of the photons is hard. The 
virtual and soft corrections to the initial-state vertex are 
identical to the muon case and Eq. (|12f> holds, with <icrp g R 
now being the leading order e + e~ — > 7r + 7r - 7 cross section, 
and the photon being emitted off the final pions only. 

We report here also formulae analogous to Eqs. (fTUf> 
and (|15|) . as we have found representations simpler than 
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Fig. 9. Virtual (a) and real (b) radiative corrections to real 
FSR emission from two-pion final state. 
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previously reported in |18| . The virtual plus soft correc- 
tions, with soft photon energy cut E% ut in the s' rest frame, 
read 



?7 V+s ( S ',£2 Ut ) = -2 
2 + Pl 



2(3, 



log(^) + 1 



\og{2w) 



(3, 

2(3, 
4Li 2 



log(t„) - 2 + log 
21og(^)log 



l-(3l 
4 

1 + Ar 



2(3, 



l + (3, 



(21) 



where w = E^/Vs*, (3, = yjl - ~Am%Js' and t, = (1 
(3,)/(\ + (3,), while, with a soft photon energy cut E, 
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in the e + e cms frame, they read 



2/3. 



x log(2w) + l + ^— log I - 
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log(£.) log 



1+Ar 



4Li s 
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(22) 



posed they contribute up to 2-3% to the Q 2 differential 
cross section, as shown in Fig. ^| However, for KLOE 
cuts used for low-angle measurements of the pion form 
factor, the additional IFSNLO corrections, not included 
in PHOKHARA 3.0, are well below 0.1%, as anticipated 
in |18|. It is shown in Fig. ^] for event selections with and 
without the track mass cut. It is worth mentioning that 
the IFSNLO corrections change also the ratio 



where w = E™ in /^/s. 

The numerical importance of the corrections from Fig-EI 
depends on the event selection. If no angular cuts are im- 



4(1 



2ml 



Iq 1 



f%\FAQ 2 



(23) 



dQ 1 



which is expected to be 1 for ISR only, when no cuts are 
applied. From Figs. 1121 and!13lit is clearly visible that one 
has to take the IFSNLO corrections into account and use 
a Monte Carlo event generator, when extracting the pion 
form factor from the ratio of pion and muon cross sections. 
This is true at both <P- and i3-factories. A similar situation 
can be expected for other hadronic final states. 



5 Vacuum polarization implementation 

PHOKHARA 4.0 includes also vacuum polarization cor- 
rections for all final hadronic states and muon production, 
with the actual implementation of Ref. [H1J ■ As can be an- 
ticipated from Fig. ^ the vacuum polarization corrections 
are calculated at different scales for different types of dia- 
grams (Q 2 , Q 2 , s, etc.). In the energy range that is impor- 
tant for measurements of the components of the hadronic 
cross section through the radiative return method, the vac- 
uum polarization corrections have a non-trivial behaviour 
(see Fig. I14|) due to the hadronic contributions. For event 
selections, for which only ISR corrections are important, 
the vacuum polarization contribution is just a multiplica- 
tive factor depending on Q 2 , so that one can correct the 
event generator results even after generation. However, for 
event selections where other contributions are also impor- 
tant, event generation with complete implementation of 
vacuum polarization might become necessary. An exam- 
ple is shown in Fig. 1151 where the vacuum polarization 
contributions do change the forward-backward asymme- 
try of the pion polar angle distribution. The effect, of up to 
5%, is mainly caused by the choice of s in the close vicin- 
ity of the <P resonance, which is nothing but the KLOE 
case. As seen in Fig. 1141 the vacuum polarization around 
the <P is much larger than at lower energies. As a result, 
the vacuum polarization correction to the ISR-FSR inter- 
ference, which contributes to the numerator of the asym- 
metry, is larger than the corresponding corrections to its 
denominator, which is dominated by ISR. Therefore, vac- 
uum polarization corrections do not drop in the ratio. This 
effect has to be taken into account, when studies of model 
dependence of FSR are performed. 



6 Summary and Conclusions 

Precision measurements of the hadronic cross section through 
the radiative return to an accuracy of about 1% require the 
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full understanding of radiative corrections and their im- 
plementation in a Monte Carlo event generator. For this 
reason a number of new features have been introduced 
into the Monte Carlo event generator PHOKHARA. 

In a first step we presented the formulae for real and 
virtual radiation, which are required to describe FSR for 
muon pair production through the radiative return in next- 
to-leading order. This includes a specific set of virtual 
corrections as well as real radiation. The Monte Carlo 
event generator PHOKHARA has been upgraded and, in 
its present form, also simulates the contribution from the 
"two-step" reaction e + e~ — » 77* (— > The techni- 

cal precision of the generator in this new mode has been 
demonstrated to be better than 10~ 4 . It is shown that 
the relative importance of this NLO correction depends 
strongly on the details of the experimental cuts. 

The program has also been extended to the case of 
pion production. In particular corrections from virtual 
plus soft ISR, in combination with hard FSR, have been 
introduced. Although the impact of this additional piece 
is small, as far as actual measurements at KLOE or the B- 
factories are concerned, it is required for a complete NLO 
treatment of the radiative return for pion-pair produc- 
tion. Another straightforward upgrade is the introduction 
of vacuum polarization in the photon propagator. As far 
as charge-symmetric cross section measurements are con- 
cerned, this effect can just be considered as a multiplica- 
tive, Q 2 -dependent change of the normalization. However, 
interference terms between amplitudes with virtual pho- 
tons of different virtualities, which arc relevant to e.g. the 
forward-backward asymmetry, do lead to modifications of 
the results at the level of several per cent. 
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